UTAP-574 



Dynamics of colliding branes and black brane production 
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We study the dynamics of colliding domain walls including self-gravity. The initial data is set 
up by applying a BPS domain wall in five-dimensional supergravity, and we evolve the system 
determining the final outcome of collisions. After a collision, a spacelike curvature singularity 
covered by a horizon is formed in the bulk, resulting in a black brane with trapped domain walls. 
This is a generic consequence of collisions, except for non-relativistic weak field cases, in which 
the walls pass through one another or multiple bounces take place without singularity formation. 
These results show that incorporating the self-gravity drastically changes a naive picture of colliding 
branes. 

PACS numbers: 98.80.Cq, 11.25.Wx 



I. INTRODUCTION 

It is known that primordial black holes and domain 
walls may have been produced in the early universe 
through the physical process of the collapse of cosmolog- 
ical density perturbations and the scries of phase tran- 
sitions during the cooling phase of universe. On the 
other hand, black holes and domain walls (also known 
as branes) also play an important role in string theory as 
fundamental constituents. In addition, according to M- 
theory, branes are of particular relevance to cosmology: 
branes arc free to move in a bulk space, and they may 
approach and collide, causing the big bang/crunch or an 
inflation on branes 

In view of the phenomenological relevance, under- 
standing how the domain walls/branes interact dynam- 
ically is an important problem, and more knowledge in 
this area could help in clarifying many issues regarding 
the early universe. In the past few years much atten- 
tion have been paid to understanding the dynamics of 
domain walls and bubbles (e.g., 0, [l] and 0]). In par- 
ticular, the interaction between black holes and domain 
walls has been the subject of study. Nevertheless, even 
more fundamental processes like collision, recoil, and re- 
conncction of branes arc less understood. 

The collision and recoil of domain walls in the cosmo- 
logical context described above was studied in , where 
a reheating mechanism via particle productions was dis- 
cussed within a toy model. In this paper, we consider the 
problem from a different perspective. The collision of do- 
main walls/branes is a violent phenomenon, and, as par- 
tially observed in our previous study, a spacctime singu- 
larity might appear through a collision. If this is the case, 
a low-energy description of colliding branes breaks down 
at some point, impling a complete loss of predictability, 
without the complete theory of quantum gravity. 

We investigate the process of collision using a BPS do- 



main wall in fivc-dimcnsional supergravity, and our main 
goal is to determine the final outcome of the kink-anti- 
kink collisions including self-gravity. As we will sec, sin- 
gularity formation is a generic consequence of collisions. 
However, the singularity is spacelike and hidden inside 
the horizon. The horizon extends in a spatially flat di- 
rection along the brane so that a black brane is produced 
through the collision. To clarify and to provide further 
examples of black brane production, we will also study 
collisions using another model of domain walls. 

II. BRANE DYNAMICS 

A. Model of virall collisions 

The system we intend to study consists of two do- 
main walls that are initially located far away from each 
other. The initial data for such a configuration is con- 
structed by superposing domain wall solutions in an ap- 
propriate manner. As a model of a domain wall, we con- 
sider a gravitationally interacting scalar field of the La- 
grangian C ~ — ^ V^^Vqi/) — V. The following solution 
of 5D Einstein equations represent a single domain wall 
{Model I), which has a spatially flat direction in three 
dimensions dxi,. 
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where cjq = 3 - 4a^L^, wi = -8aL^, UJ2 = -3 - AL"^. 
There are three unfixed parameters, i.e., wall thickness 
5, amplitude L of scalar field, and the position tq of the 
kink's core. We will hereafter take = 1. We will 
call this domain wall solution the kink solution (for a < 
0). The anti-kink solution is defined by the reflecting 
r-coordinate in the above solution. 
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(a) 5k = 5a = 1, Lk = La =0.18 



(b) 5k = 1 and 5a 



(c) Lk = 0.18 and La = 0.22 



FIG. 1: Energy density pont = const, surfaces: (a) collision of two identical walls at the center, (b) collision of two walls with 
different thickness, and (c) two walls with different amplitude. The sharp peaks represent the domain walls, and the arrows 
show the directions of a wall's velocity. For (b) and (c), all the unspecified parameters are the same as those in (a). 
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FIG. 2: R and R"'^'^''' Rated at t =const. surface, corresponding to the simulations in Fig. [T] 



In the limit of r ^ ±oo, the scalar field asymptotes 
constants, and the scalar potential plays the role of the 
cosmological constant A = k^T^ in the limit 



A = --^[l + « sign(r)J^ 



■+±oo 



(2) 



The domain wall for |a| < 1 gives a warp factor decreas- 
ing for both infinities of an extra dimension, and the 
cases of |a| = 1 become the wall solutions interpolating 
between AdS and flat Minkowski vacua. For \a\ > 1, the 
warp factor decreases in one direction, and increases in 
the other. This domain wall can be embedded into the 
five-dimensional supergravity coupled with hypermulti- 
plets as an exact BPS domain wall (B]. After integrat- 
ing out irrelevant fields with canonical normalization, the 
above solution is found to be identical to the exact BPS 
solution in 0]. 

We shall restrict our analysis to collisions along a r- 
direction, preserving the symmetry along the homoge- 
neous afs-directions. Even with this simplification, such 
a setup is of relevance in a number of physical situations. 
The initial data for such a collision can be obtained as 
follows. First of all, we introduce a new coordinate z by 
z = J dre~^ and work on the conformal gauge. 



(3) 



Then the above single static wall is boosted along the fifth 



direction z, and we obtain a wall moving with constant 
velocity (3 

To discuss collisions of two moving domain walls, we 
set a kink solution at z = — zq and an anti-kink solution 
at z = zq, which are separated by a large distance and 
approaching each other with the same (or different) speed 
(3. Such superposition and matching of the metric and 
scalar field at the center is possible for \a\ = 1, and suf- 
ficiently smooth initial data that satisfies the constraint 
equations at the initial time can be obtained, as long 
as the spatial separation between the two walls is much 
larger than the thickness of walls. Therefore, we take 
\a\ = 1 throughout this paper. Obviously, we set A = B 
and velocity A = B at the outset, and the initial values of 
(j) and A are given by the above construction. During the 
evolutions, the Neumann boundary conditions are im- 
posed at the outer boundaries. The asymptotics of the 
scalar field is oc (7|z|/y67|A| -f and the met- 

ric behaves oc {j\z\/ •\/6/|A| + 1)^-'- as |z| oo, where 
7 = 1/a/1 — is the Lorentz factor. The kink and anti- 
kink solutions are characterized by their own width and 
amplitude, {Sk,Lx) and ((5^,L^), respectively. There- 
fore, we have three types of unfixed parameters for the 
initial setup; S, L and /3. 

Using a fourth-order accurate finite difference code, we 
have solved the system numerically and evaluated the 
constraints at each time step for various families of initial 
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FIG. 3: (a) and (b) show the numerical results in the null coordinates, corresponding to Fig. 1(a) and 1(b) respectively. 
Dotted lines with arrows describe the orbit of walls. Thick lines are u, v-ax.es, and dotted straight lines show u,v =const., to 
which the spacelike singularity asymptotes. The position of a wall is defined by its maximum energy density, (c) Time variation 
of the (kink) wall's speed. Irrespective of the initial speed, the final speed after the collision goes to unity. 



data. The overall picture does not depend on a specific 
choice of the parameters. Some examples of numerical 
results are reported in Fig.[l] in which the time evolutions 
of energy density 



V 



(4) 



in {t, z}-flamc are shown. In all these cases, the two 
walls with initial velocity (3q = 0.4 collide at z = and 
t K, 31. Fig. 1(a) describes the symmetric collision of two 
identical walls. In this case, the walls pass through one 
another so that the initial kink solution at z < goes 
to z > 0. (The kink and anti-kink solutions are distin- 
guished by their relative field values at the center and 
infinity z ~ ±oo.) The energy density of wall, i.e. the 
wall's tension, increases during this process. This would 
be caused by the fact that the induced universe on the 
walls are contracting during the process, with B < 0. 
After the collision, a sharp peak of density appears at 
the collision point z — Q, and it implies an emergence of 
singularity. In fact, the curvature diverges rapidly at the 
point, whereas the curvature on the wall remains finite 
and small at the moment (Fig. Here our criterion 
of curvature singularity is that the Kretschmann scalar 
exceeds R'''""^Rabcd > 10^ At the time t « 42 of singu- 
larity formation, the energy density localized at z = is 
1.2 times bigger than those on the walls, and a portion of 
energy is stored in this small region, which will be inside 
an event horizon, as we see below. 

This basic picture of collision holds for other cases. For 
the asymmetric collisions, such as two walls with differ- 
ent width, amplitude, and/or speeds, the emergence of 
singularity is still a generic feature. Figs. 1(b) and 1(c) 
show examples of collisions in which different thickness 
or amplitude of scalar field are taken for the two walls, 
without changing other parameters. Among these cases, 
Fig. 1(b) shows that one of the walls recoils at the col- 



lision, due to the larger momentum of one of the walls: 
the initial kink solution at z = —30 in Fig. 1(b) goes to 



z > after the collision, while the anti-kink at z = +30 
bounds back. Interestingly, in these asymmetric colli- 
sions, the curvature singularities appear off the collision 
point. For Figs. 1(b) and 1(c) they are at z = —2.4 and 
z = —1.2, respectively. 

For the wide range of initial parameters, the emergence 
of singularity is the generic consequence. However, as 
expected and discussed below, the singularity does not 
appear for L 1 and/or /3o ^ 1 for fixed 6. In such 
"non-relativistic" cases, the two walls just pass through, 
and the final configurations of fields are well described 
by the boosted walls, as we applied for the initial config- 
urations. 

B. Horizon formation 

The next task at hand is to confirm the nature of sin- 
gularity. In numerical investigations of singularity for- 
mation and global structure of the spacctime, null co- 
ordinates are useful to prevent the singularity from cor- 
rupting the rest of the spacetime. In these coordinates, 
horizons are not particularly special and we can follow 
the collision all the way to the singularity even when a 
horizon appears through a collision. We evolve the col- 
liding walls in the double-null coordinates (e.g., [1, 0? 



ds'' = -2e'^^dudv + 



„2B 



(5) 



where ^/2u = {t — z), \/2v = {t + z). In this gauge, 
the Einstein equations and the dynamical equation for a 
scalar field are split into three dynamical and two con- 
straint equations. 

Let us first focus on the symmetric collision in 



Fig. 1(a) The corresponding evolution in the null coor- 



dinates is described in Fig.[3l It shows that the curvature 
singularity is spacelikc, approaching u or v = const, lines 
at late times, which corresponds to the event horizon. 
This result is very generic, and wc have observed similar 
results for the wide range of initial data (velocity, etc.). 



4 



Minkowski 



AdS 





Kink 



Antikink 



FIG. 4: Schematic conformal diagrams for a single domain 
wall that asymptotes to AdS (Minkowski) as z oo {z 0) 
(Left) and for colliding walls, producing a black brane (Right). 



This system has homogeneous 3-spatial directions, and so 
the horizon also extends in these directions. This means 
that a black brane is produced by the collision of walls, 
so that this collision provides the dynamical mechanism 
of generating black branes in higher dimensions. 

Another interesting feature is that after the collision 
the walls are trapped around the surface of the horizon. 
The final speeds of walls asymptote to the speed of light 
[Fig. 3(c)| , irrespective of the initial velocity. The bulk 
outside the two walls is not exactly AdS, but asymptotes 
to it. Because of this behavior, the walls are pulled out- 
side, accelerating in the directions. A schematic picture 
of a conformal diagram is given in Fig. [4] 

Another example corresponding to to the asymmetric 

Even in this 



collision in Fig. 1(b) is shown in Fig. 3(b) 



asymmetric collision, the event horizon forms from the 
point where the spacelikc singularity appears. An inter- 
esting difference is that the kink wall escapes from the 
horizon, and only the antikink wall is trapped nearby. 

C. Model II 

The initial data discussed so far is based on the sin- 
gle BPS domain wall. There is another simple model in 
which collisions of walls can be tested. It is the model 
used in the previous work 01, 



U 
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where W ~ (f'^) vfl'^^ I' ^^'^ basic property 

of the wall is quite similar to the wall in the previous 
sections; the bulk in r > asymptotes to the Minkowski 
spacetimes, while the spacetime in r < asymptotes to 
the AdS, recovering ([2]). This single domain wall solu- 
tion is found simply by extending the four-dimensional 
solution in M- 



It is interesting to study various aspects of the wall 
collisions in this model and compare them with the pre- 
vious model. We have performed many simulations and 
confirmed that all phenomena observed in the previous 
sections, such as the singularity and horizon formation 
hold with qualitatively similar behaviors. A basic ex- 
ception is that in this model the walls bounce back after 
the collision, contrary to the case in Model I. Thus the 



causal structures of such a collision look like Fig. 3(a) 



but kink and anti-kink profiles are exchanged after a col- 
lision. This difference comes from the nonlinear interac- 
tion through the dynamics of collision. 

Such details of model dependence become more signifi- 
cant for weak field cases in which no singularity appears. 
In Fig.O we compare the difference of the two models by 
showing the orbits of a wall after the symmetric collisions 
for various values of L. For Model I, the spacelikc sin- 
gularity appears for L > 0.1, and the walls asymptote to 
the null lines, as discussed above. On the other hand, for 
L < 0.1, the velocity of the wall becomes timelike with 
constant speed after the collisions, and no singularity ap- 
pears. In fact, the final configurations of scalar fields 
are well approximated by superposing boosted walls, so 
that the two walls just pass through one another in these 
cases. Note that if the initial velocity is increased the 
horizon appears even for smaller L. 

For Model II, multiple collisions take place for L <C 1 
(Fig. [5]). For L = 0.01, the collision takes place two 
times, and then the wall bounces back with constant ve- 
locity. This behavior is compatible with and typical in 
the non-gravitating system of the previous study ■ As 
L increases, the two walls gravitate toward one another 
and multiple bounces take place (e.g., L = 0.045). The 
marginal value of L is L = 0.05 in Fig. [51 in which the 
gravitational attractive force and the repulsive force due 
to outer AdS region are in balance. Therefore, a quasi- 
static configuration of two walls is realized after the col- 
lision. 

III. SUMMARY AND DISCUSSION 

We have considered a system of colliding domain walls, 
whose initial data is set up based on a single BPS domain 
wall, and found that the role of gravity at the collision 
is significant in that it can drastically change the pic- 
ture of "silent" collisions without self-gravity, as observed 
also in colliding (impulsive) gravitational waves [ll| . The 
main result of our study is that horizon formation is a 
generic phenomena in the collision of walls. In the non- 
relativistic cases, such as L <C 1 and/or /? <C 1, "silent" 
collisions without singularity and horizon are realized, 
but such cases are very limited and unlikely in an early 
universe of braneworlds. 

The local interaction and dynamics at the collision de- 
pends on which model we are looking at. For Model I, 
the domain walls can pass through one another, while 
in Model II they bounce at the collision and go back. 
However, the basic feature of horizon formation does not 
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FIG. 5: The orbits of a wall after the symmetric collisions at 
t 60 (/3o = 0.2 and 5 = 1) are described for the two models 
(Model I, II), to show the dependence on L. Solid straight 
(black) lines represent u =const. lines along the spacelike 
singularities for L = 0.1, 0.2. 



change in these models, and we have concluded that the 
horizon formation (and singularity behind it) is a generic 
consequence of kink-anti-kink collisions. The horizon has 
three homogeneous spatial directions, so that a black 
three-brane is produced by such a collision. The bulk 
outside the two walls is approximated by the AdS, and 
then we look at the created black hole as it sits on the 
AdS. In the pure AdS, a possible black hole is a topo- 
logical black hole, which has a flat 3-dimensional hyper- 
surface with vanishing curvature. The end state of the 



present scenario would be this type of black hole (Fig. [4]), 
although the field fills the bulk outside the horizon. Here, 
a further interesting possibility comes from the fact that 
a spatially homogeneous horizon suffers from Gregory- 
Laflammc instability in general. The end state of this in- 
stability has not been clearly understood so far, and the 
horizon may break up, resulting in multiple black holes 
that are stuck on the walls Q. Furthermore, during this 
process, a good deal of energy will be radiated away by 
gravitational waves, and they remain as primordial grav- 
itational wave backgrounds. Thus this possibility pro- 
vides a new way of producing primordial black holes and 
gravitational waves in an early universe with higher di- 
mensional bulk filled by walls/branes. This scenario is 
analogous to the bubble collisions in a four-dimcnsional 
inflationary universe [l^ . 

There are several interesting directions, which may be 
pursued on the basis of these results and/or by relaxing 
several conditions. One of such issue is a study of the 
effects and roles of other fields in supergravity, a subject 
sets aside in our analysis. Other fields contained in the 
hypcrmultiplcts will be excited (or thcrmalizcd) during 
the collisions, and analyzing them should provide many 
cosmological insights into braneworld cosmologies [1C| . 
Furthermore, there is one most important question left 
to answer: To what extent is the production of black 
holes/branes generic in a class of more generic theory 
and context, such as collisions of different types of walls 
with an arbitrary incident angle. We foresee no major 
obstacle in anticipating that horizon formation would be 
suppressed. 
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